Study of the superconducting state in the Cmmm phase of GeH4 compound 
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The superconducting state, that gets induced in the GeH4 compound under the pressure at 20 
GPa, is characterized by the high critical temperature Tc (li*) £ (54.5,31.2) K, where fj,* is the 
Coulomb pseudopotential, and /x* G (0.1, 0.3). Other thermodynamic parameters for the considered 
range of the Coulomb pseudopotential significantly differ from the predictions of the BCS model. 
In particular: (i) the ratio of the energy gap to the critical temperature (Ra = 2A (0) /kBTc) 
changes from 4.07 to 3.88; (ii) the value of the parameter Rc = AC (Tc) jC^ (Tc) decreases from 
1.85 to 1.69, where AC denotes the specific heat jump, and C^ is the specific heat of the normal 
state; (iii) the ratio Rh = TcC^ (Tc) /He (0) increases from 0.147 to 0.155; He represents the 
thermodynamic critical field. 
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In the elements and the chemical compounds under the 
influence of high pressure (p), the electron-phonon inter- 
action can be so strong that the superconducting state 
is characterized by a high value of the critical temper- 
ature. Of particular note are the rcsuhs obtained for 
hthium and calcium, where the maximum value of the 
critical temperature is equal to: [TclpLao 2GPa = 14 K 

and [Tclp^ieiGPa = 25 K Q, [2]. The other thermody- 
namic parameters of the superconducting state in Li and 
Ca, as the result of the existence of the retardation and 
strong-coupling effects, cannot be properly described in 
the framework of the BCS model [|, 0, j|, 0, 0- 

Probably, the highest value of the critical temperature 
can be observed in metallic hydrogen [§] . Unfortunately, 
also at an extremely high pressure. In particular, it has 
been predicted that hydrogen enters the metallic state 
at 400 GPa [§]. It exists then in the molecular state, in 
which the superconducting state with the critical tem- 
perature of 240 K has been predicted [l^l- Note, that 
the superconducting state in the molecular hydrogen can 
be strongly anisotropic in some ranges of the pressure. 
In the simplest case, the three-band formalism should be 
used for the description of its properties [l^l , [HI ■ 

Above 500 GPa, molecular hydrogen dissociates into a 
single atomic phase [l^, 13 1. In the examined case, the 



highest value of Tc is predicted for the pressure at 2000 
GPa, where the critical temperature reaches 600 K [l3| . 
For the considered pressure, the low temperature energy 



gap at the Fermi surface (2A (0)) is so large that the 
dimensionless ratio 2A (0) /kgTc equals approximately 
6 [m . The comparable value of 2 A (0) is only observed 
in the group of the high-temperature superconductors 
(cuprates) [ll, [l3|. 

In order to reduce the pressure of metallization in hy- 
drogen, a method of the chemical precompression has 
been proposed [l^l- Physically, this idea can be realized 
in the group of the hydrogenated compounds of the fol- 
lowing_Upe: SiH4 {^2)2: ^^2^6, B,H6 and PtH [13, [13], 

HI, M- 

Typically, the properties of the superconducting state 
for pressures that excess the value of 100 GPa have been 
considered in the literature. Hence, it is difficult to verify 
the theoretical results in the experimental way. 

In the paper, we have determined the basic thermo- 
dynamic properties of the superconducting state in the 
Cmmm phase of the GeH4 compound for the relatively 
low value of the pressure {p = 20 GPa). In the examined 
case, the critical temperature, calculated on the basis of 
the modified McMillan formula, equals ~ 47 K [H, [H], 
[25| . Let us note that the low value of p will allow in the 
future an easy comparison of the results obtained in this 
study with the experimental data. 

Numerical calculations have been conducted in the 
framework of the Eliashberg formalism. On the i mag - 
inary axis the Eliashberg equations take the form j2a |. 
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where </>„=(/) (icLi„) denotes the order parameter func- 
tion, and Z„ = Z is the wave function renormal- 
ization factor; aj„ stands for the Matsubara frequency: 
LOrn = (tt//?) (277Z — 1), where l3 = {kBT)~^; ks is the 
Bohzmann constant. The order parameter is defined by 
the formula: A„ = 0„/Z„. 

The pairing kernel for the electron-phonon interaction 
can be represented by: A (z) = 2 /g^"""' dft F (17), 

where the Eliashberg function {a^F (fi)) has been de- 
termined in the paper (23| . The value of the maximum 



J 



phonon frequency (rJmax) is equal to 515.69 meV. 

The symbol fi* represents the Coulomb pscudopoten- 
tial, which models the depairing electron-electron inter- 
action. The quantity 9 denotes the Heaviside function, 
and cUc is the cut-off frequency: uJc = Sfimax- 

The Eliashberg equations have been solved with the 
help of the numerical methods discussed and used in [28j . 
The convergence of the solutions has been obtained in the 
temperature range from Tq = 3 K to Tc; M = 1100 has 
been assumed. 

The form of the functions (f> and Z on the real axis (uj) 
has been determined by using the Eliashberg equations 
in the mixed representation [29|: 
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where the symbols N [uj) and / {uj) represent the Bose- 
Einstein and Fermi-Dirac functions, respectively. 

A way of solving the equations ^ and ^ has been 
discussed in [s^l ■ 

Fig.[l]represents the dependence of the critical temper- 
ature on the Coulomb pseudopotential {^* G (0.1,0.3)). 
The obtained results prove that Tq is high in the whole 
range of /i* which has been considered. In particular: 
Tc e (54.5,31.2) K. Additionally in Fig.[Tl we have plot- 
ted the courses of the dependence of Tc on ^* obtained 
with the help of the McMillan and Allen-Dynes expres- 
sion [i^, [2^. It can be easily seen that the analytical 
formulas significantly underestimate the critical temper- 
ature, especially for the high values of the Coulomb pseu- 



dopotential. 

In Fig.O (A)-(C), we have plotted the order param- 
eter on the imaginary axis for the selected values of 
the temperature and the Coulomb pseudopotential. It 
has been found that with the increase in the num- 
ber m, the function A™ strongly decreases, and then 
becomes saturated. The influence of the temperature 
and the Coulomb pseudopotential on the order param- 
eter value can be the most conveniently traced on the 
example of the course of the function A,„=i (T, /x*). 
It can be parameterized with the help of the for- 



mula: A„=i (T, n*) 



A„=i(ro,M*)\/i 



where A„=i (Tq, /x*) = 58.4 (^*)-44.58/i^ 4-13.286 mcV 
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FIG. 2: (A)-(C) The form of the order parameter on the 
imaginary axis for the selected values of the temperature and 
the Coulomb pseudopotential. (D) The order parameter for 
m = 1 as the function of the temperature for the selected 
values of the Coulomb pseudopotential. 



and /3 = 3.3. 

It should be noted that the quantity 2A„j=i (T, /i*) 
with a good approximation reproduces the dependence of 
the energy gap at the Fermi surface on the temperature 
and the parameter /i*. The obtained results have been 
shown in Fig.[2](D). 

The second solution of the Eliashberg equations on the 



FIG. 3: (A)-(C) The form of the wave function renormal- 
ization factor on the imaginary axis for the selected values 
of the temperature and the Coulomb pseudopotential. (D) 
The renormalization factor for m = 1 as the function of the 
temperature for the selected values of the Coulomb pseudopo- 
tential. 



imaginary axis has been presented in Fig. [3] (A)- (C). Just 
as it was in the case of the order parameter, the value 
of Zm rapidly decreases with the increasing number m. 
On the other hand, the effect of the temperature and the 
Coulomb pseudopotential on the wave function renormal- 
ization factor is meager. 

It is worth noting that for T = Tc, the value of Zm=i 
is independent of the parameter /i*. In the considered 
case, it can be calculated analytically using the formula: 
Zm=i = 1 -t- A, where the symbol A is the electron- 
phonon coupling constant: A = 2 J^"""" (f^) F (il) 
For the GeH4 compound A is equal to 1.0324, thus: 
[Zm=i]rp^rp^ = 2.0324. The identical result has been ob- 
tained by using the numerical procedures, indicating a 
high accuracy of the computation carried out. 

From the physical point of view, the wave function 
renormalization factor for m = 1 with a very good ap- 
proximation reproduces the value of the electron effective 
mass: m* ~ Zm^irue, where the symbol rrie denotes the 
electron band mass. The results presented in Fig.|3](D) 
prove that the electron effective mass in GeH4 is high 
in the entire range of existence of the superconducting 
state. 

The solutions of the Eliashberg equations on the imag- 
inary axis have been used in the calculations of the free 
energy difference between the superconducting and the 
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The symbols and denote the wave function renor- 
malization factor for the superconducting state {S) and 
the normal state {N), respectively. 
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FIG. 4: The free energy difference (at the bottom) and the 
thermodynamic critical field (at the top) as the function of 
the temperature for the selected values of the Coulomb pseu- 
dopotential. 
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FIG. 5: The specific heat of the superconducting state and 
the specific heat of the normal state as the functions of the 
temperature for the selected values of the Coulomb pseudopo- 
tential. 

In Fig. 131 below the value, there has been pre- 
sented the dependence of the free energy difference 
on the temperature for the selected values of the 
Coulomb pseudopotential. It can be easily seen that 
with the increase of ji* the absolute values of the 
free energy strongly decrease. In particular, the ratio 
[AF (To)]^.^o.i / [AF (To)]^.^o.3 is equal to 3.20. 



The thermodynamic critical field has been calculated 
by the means of the formula: 



Hr 



v/-87r [AF/p(0)]. 



(6) 



Above the value in Fig. |4] obtained course have been 
plotted. Also in this case, the depairing electron correla- 
tions significantly reduce the values of the function un- 
der consideration: [He (0)]^.^o ^ / [Hq (0)]^,^(, 3 = 1.77, 
where the designation He (0) = He (To) has been intro- 
duced. 

The specific heat difference (AC = — C^) has been 
estimated on the basis of the formula: 



AC (T) 



\ <P [AF/p(0)] 

/3 



(7) 



kBp{Q) P d{kBTY 

On the other hand, the specific heat of the normal state 
may be calculated by using the formula: kspio) ~ ^' 
where the Sommerfeld constant is derived from: 7 = 
|7r2(l + A). 

The results obtained for and have been pre- 
sented in Fig.[S] It can be seen that at higher tem- 
peratures the specific heat of the superconducting state 
far exceeds the value of , which results in the jump 
at the critical temperature. In addition, we can point 
out that the specific heat jump strongly depends on 
the assumed value of the Coulomb pseudopotential: 
[AC (Tc)]^.^o.i / [AC (Tc)]^.^o.3 = 1.92. 

The considered thermodynamic functions allow us to 
calculate the values of two dimensionless ratios. 



Rh 



_ TcC^ [Tc) 



HhiO) 



and 



AC(rc) 



C^ (Tc) ■ 



(8) 



In the framework of the BCS theory, their values 
are universal and respectively equal to: 0.168 and 1.43 
0i 0- For Gell4 compound, it has been obtained: 
Rh e (0.147,0.155) and Rc e (1.85,1.69), when fj.* e 
(0.1,0.3). On this basis, it is clear that the thermody- 
namic parameters of the superconducting state in Gell4 
significantly deviate from the expectations of the BCS 
model. 

The physical value of the energy gap at the Fermi level 
should be determined on the basis of the expression: 



A (T) = Re [A (w = A (T))] 



(9) 



The courses of the order parameter function on the 
real axis have been obtained by solving the Eliashberg 
equations in the mixed representation. The results have 
been presented in Fig.|6l 
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FIG. 6: The real and the imaginary part of the order parameter on the real axis for the selected values of the temperature and 
the Coulomb pseudopotential. Additionally, the rescaled Eliashberg function has been plotted. 



It can be seen that the order parameter takes com- 
plex values. For the lower frequencies, in the area of 
the applicability of the equation © , it can be conducted 
that: Re [A (oj)] ^ Im [A (w)]. From the physical point 
of view, this result indicates the existence of the weak 
damping effects [s^]- In addition, we can note that the 
complex behavior of A (w), in the frequency range from 
to about 300 meV, is explicitly related to the shape 
of the Eliashberg function. Above 300 meV the order 
parameter becomes saturated. 

On the basis of Fig. [HI we can also trace the influence 
of the temperature on the course of A (w). It has been 
found that together with the increase of T, the outline 
of the order parameter function becomes more regular, 
regardless of the value assumed for the Coulomb pseu- 
dopotential. 

The presented results allowed us to calculate the value 
of the energy gap at the Fermi level for the lowest consid- 
ered temperature (2 A (0)). Then, the dimensionlcss ratio 
has been estimated: Ra = 2A (0) /kBTc- In the case of 
the GeH4 compound, the following result has been ob- 
tained: Ra G (4.07,3.88), for G (0.1,0.3). Let us 
notice that the BCS theory predicts that: Ra = 3.53. 
Therefore, the parameter Ra for GeH4 clearly exceeds 
the value estimated in the framework of the BCS model. 

To summarize, we have calculated the values of the 
basic thermodynamic parameters of the superconducting 



state in the Gell4 compound. The calculations have been 
performed for a wide range of the Coulomb pseudopoten- 
tial: ij* G (0.1,0.3). It has been found that even for the 
large fj.* , the critical temperature reaches a relatively high 
value (see table [T|. The values of the dimensionlcss ther- 
modynamic ratios (i?Ai Rc ^-nd Rh) significantly differ 
from the predictions of the BCS theory. However, with 
the increasing value of the Coulomb pseudopotential, the 
difference between the Eliashberg and BCS model fades 
(table H. 

It should be also noted that the value of the spe- 
cific heat jump at the critical temperature and the low- 
temperature value of the thermodynamic critical field de- 
creases strongly with the increasing of fj.* , which has been 
explicitly presented in table HI 
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TABLE I: The values of the thermodynamic parameters of the superconducting state in the GeH4 compound for p = 20 GPa. 







Tc (K) 


-Ra, -Rc, Rh 


AC{Tc) /kBp{Q) (meV) 


i/c(0)/Vp(0) (meV) 




= 0.1 


54.5 


4.07, 1.85, 0.147 


116.31 


44.80 




= 0.2 


39.9 


3.95, 1.79, 0.151 


82.45 


32.30 




= 0.3 


31.2 


3.88, 1.69, 0.155 


60.72 


25.03 
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